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SGR J1550-5418 exhibited three active bursting episodes in 2008 October and in 2009 January 
and March, emitting hundreds of typical Soft Gamma Repeater (SGR) bursts in soft gamma-rays. 
The second episode was especially intense, and our untriggered burst search on Fermi/GUM data 
(8—1000 keV) revealed ~450 bursts emitted over 24 hours during the peak of this activity. Using 
the GBM data, we identified a ~150-s-long enhanced persistent emission during 2009 January 22 
that exhibited intriguing timing and spectral properties: (i) clear pulsations up to ~110 keV at the 
spin period of the neutron star, (ii) an additional (to a power-law) blackbody component required 
for the enhanced emission spectra with kT ~ 17 keV, (iii) pulsed fraction that is strongly energy 
dependent and highest in the 50—74 keV energy band. A total isotropic-equivalent energy emitted 
during this enhanced emission is estimated to be 2.9 x 10 40 (D/5 kpc) 2 erg. The estimated area of 
the blackbody emitting region of « 0.046(D/5 kpc) 2 km 2 is the smallest "hot spot" ever measured 
for a magnetar and most likely corresponds to the size of magnetically-confined plasma near the 
neutron star surface. 



I. INTRODUCTION 

IE 1547.0-5408 was observed with the X-ray Multi- 
Mirror Mission (XMM-Newton) in 2004 as a magne- 
tar candidate, selected for its galactic plane location 
and its relatively soft magnetar-like spectrum as seen 
with the Advanced Satellite for Cosmology and Astro- 
physics (ASCA) during their Galactic plane survey 
Although no period was detected in the original and 
follow-up XMM observations, its positional coinci- 
dence with an extended galactic radio source G327.24- 
0.13 (possibly a supernova remnant) also suggested 
that IE 1547.0-5408 is a magnetar 0. The subse- 
quent discovery in radio observations of a spin period 
of 2.07 s and a period derivative of 2.3 x 10~ n s s _1 
led to an estimated dipole surface field of B ~ 2.2 x 
10 14 G and confirmed the source's magnetar nature; 
the source was also renamed as PSR J1550-5418 
Its period makes IE 1547.0-5408 the fastest rotating 
magnetar; the source is also one of the only two that 
emit in radio wavelengths [3, Q ■ The distance of the 
source has been estimated to be ^4—9 kpc using var- 
ious methods 0, 9, including the most recent esti- 
mate of 4—5 kpc [6[. Throughout this paper we use 
£>5 = D/5 kpc as the source distance measure. 

On 2008 October 3, IE 1547.0-5408 entered an 
episode of X-ray activity, emitting several typical 
SGR-like bursts over the next 7 days. During this 



period, 22 short duration bursts were observed with 
the Gamma-ray Burst Monitor (GBM) on board the 
Fermi Gamma-ray Space Telescope 0]. 

The source entered a second period of extremely 
high X-ray burst activity on 2009 January 22 [8] . Dur- 
ing the first 24 hours of this bursting episode, the 
Fermi/GBM triggered on the source 41 times: the 
number of triggers was limited only by the instru- 
ment's capability and did not reflect the actual num- 
ber of bursts emitted by the source. In fact, our on- 
ground search for untriggered events revealed a total 
of ^450 bursts during this 24 hour period: an un- 
usually high burst frequency from a single source [9]. 
Based on this SGR-like behaviour, we renamed the 
source as SGR J1550-5418 0. 

We discovered an enhancement of the persistent 
emission lasting ^150 s by examining the data from 
the first GBM trigger on January 22. Closer inspec- 
tion of this enhancement in various energy ranges re- 
vealed periodic oscillations with a period consistent 
with the spin period of SGR J1550-5418. Here we 
present a detailed temporal and spectral analysis of 
this enhanced emission. In i fiTl we briefly describe 
our observations and the GBM instrument and data 
types. We present our temporal analysis results in 
mill and our spectral studies in mV\ Finally we dis- 
cuss the physical implications of our discovery in i fVl 
Complete descriptions of our analysis and discussion 
are presented in Kaneko et al. (2010) (Tlj . 
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II. INSTRUMENTATION AND DATA 

In trigger mode, GBM provides three types of data 
for each of 12 Nal and 2 BGO detectors; CTIME 
Burst .CSPEC Burst, and Time Tagged Event (TTE) 
data [ll]. The CTIME Burst data have a time resolu- 
tion of 64 ms with rather coarse spectral information 
(8 energy channels). The CSPEC Burst data pro- 
vide high-resolution spectra (128 energy channels) col- 
lected every 1.024 s. Both CTIME Burst and CSPEC 
Burst accumulate data for ^600 s after a trigger. The 
TTE data provide time-tagged photon event lists for 
an accumulation time of 330 s, starting 30 s prior to 
the trigger time; this data type provides a superior 
temporal resolution down to 2/zs at the same spectral 
resolution as the CSPEC Burst data. 

The first GBM trigger at the onset of the second ac- 
tive episode from SGR J1550-5418 was on 2009 Jan- 
uary 22 at 00:53:52.17 UT (= T , GBM trigger num- 
ber 090122037). In the 600 s of the trigger readout 
we detected many individual short bursts using our 
on-ground untriggered burst search algorithm. To ac- 
cept an event as an untriggered burst, we required 
excess count rates of at least 5.5a and 4.5tr in the 
first and second brightest detectors, respectively, in 
the 10-300 keV energy range. We used CTIME data 
in both continuous (256 ms time resolution) and Burst 
mode (64 ms resolution). Subsequently, we inspected 
energy-resolved burst morphology and compared each 
detector zenith angle to the source for all 12 detec- 
tors, to determine whether the events originated from 
SGR J1550-5418. In total we identified about a dozen 
very bright bursts and over 40 less intense bursts 
within 600 s after To (see Figure [T|). During the same 
trigger readout we also discovered an enhancement in 
the underlying persistent emission starting at approx- 
imately To+70 s and lasting for ^150 s (see inset of 
Figure [T]) . 

In the analysis presented here, we have exclusively 
used data from NalO, which initially had the small- 
est detector zenith angle to the source (15°) and to 
which the source was visible through most of the en- 
hanced emission. We also note that we checked the 
LAT data (20 MeV-300 GeV) of the entire day for as- 
sociated high-energy gamma-ray emission, but found 
no evidence of high-energy photons originating from 
the direction of SGR J1550-5418. 



III. TEMPORAL PROPERTIES 
A. Timing Analysis 

During our search for untriggered events in the 
first trigger interval of 2009 January 22 from 
SGR J1550-5418, we found strong apparent periodic 
modulations in the enhanced emission period from 
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FIG. 1: GBM lightcurve of SGR J1550-5418 in 
12—293 keV. An enlarged view of the pulsed, enhanced 
emission is shown in the inset. The dashed line indicates 
the background level. 



Xo+130 to 160 s in the 50-102 keV data of detector 
NalO (see panel (c) of Figure [2]). This is the first time 
to our knowledge that pulsations unrelated to a giant 
flare from a magnetar were clearly seen in the persis- 
tent emission of an SGR, in energies up to 100 keV. 
To search for a coherent pulse period, we performed 
a timing analysis over the entire enhancement inter- 
val. We first eliminated the times of all the bursts 
found via our untriggered burst search and converted 
the remaining burst-free intervals to the solar system 
barycenter. We then generated a Lomb-Scargle peri- 
odogram [l]| [l4[ over a range of periods from 0.1 s 
to 10 s using CTIME Burst data in the 50-102 keV 
band. We found a very significant signal with a 
Lomb power of 72.6 (chance occurrence probability, 
P c ~ 10~ 16 ) at a period of 2.0699 ± 0.0024 s, which 
is consistent with the spin period of SGR J1550-5418. 
Further, to confirm our detection, we also employed 
the 7? m test (with m = 2) 15] on the burst-free and 
barycentered TTE data. We find a coherent signal 
(with a Z 2 m= 2 power of 266, P c ~ 10 -23 ) at the same 
period. Our spin period measurement is consistent 
with the one found for SGR J1550-5418 using con- 
temporaneous X-ray data (Swift/XKT) [l7| and 
radio data [Hj]. Therefore, we clearly confirm with 
the detection of these hard X-ray pulsations that the 
enhanced persistent emission seen in the inset of Fig- 
ure [1] originates from SGR J1550-5418. 

Next we searched in the enhanced persistent emis- 
sion for evolution in the intensity of the pulsations 
using a sliding boxcar technique. We found that the 
pulsed signal peaks over a 90 s interval, from To+120 
to 210 s, which encompasses the peak of the enhance- 
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FIG. 2: Lightcurve of SGR J1550-5418 in various energy 
ranges; (a) 12-27 keV, (b) 27-50 keV, (c) 50-102 keV, 
and (d) 102—293 keV. The pulsations are most prominent 
between 50—102 keV (panel c, starting at ~130 s). The 
bursts have not been removed here from the data. The 
dashed lines in panels (a) to (c) indicate the times of the 
pulse maxima. 



ment. Finally, we searched for any other intervals ex- 
hibiting pulsed emission in the burst-free continuous 
CTIME data of 2009 January 22 and during the four 
subsequent days, using a sliding boxcar of 120 s with 
10 s steps. We did not find any additional statistically- 
significant pulsed emission. For the entire search and 
for all the timing analysis reported here, we used more 
precise spin ephemeris obtained by contemporaneous 
Swift/XRT Chandra, XMM-Newton and Suzaku ob- 
servations [19j j . 



B. Pulse Profiles 

To investigate the evolution of the pulse profiles 
with energy, we folded the burst-free TTE data span- 
ning 120 s (from T o +100 s to T +220 s, which in- 
cludes the strongest pulsation period as found above) 
with the spin ephemeris of SGR J1550-5418. We es- 
timated the background level using the data segment 
between T to Xo+60 s. Figure [3] shows the source 
pulse profiles during the enhanced emission interval 



FIG. 3: Pulse profiles of SGR J1550-5418 in equal loga- 
rithmic energy intervals; (a) 10—14 keV, (b) 14—22 keV, 
(c) 22-33 keV, (d) 33-50 keV, (e) 50-74 keV, and (f) 
74—110 keV. Two cycles are plotted for clarity. The ver- 
tical dotted and dashed lines are explained in EjTlI 13 



in six energy bands that have the same logarithmic 
width. The pulse profiles above 110 keV are consis- 
tent with random fluctuations, and thus not shown. 

Figure [3] indicates that the SGR J1550-5418 pulse 
profiles in the three lowest energy bands are most 
likely complex (multi-peaked). While the two low- 
est energy band profiles are dominated by the struc- 
ture around phase 0.7-0.8 (indicated by the dotted 
lines in Figure [3]), in the 14—22 keV band we see the 
emergence of another structure around phase 0.0 (in- 
dicated by the dashed lines in Figure [3J . This pulse 
becomes equally prominent in the 22—33 keV range 
and then dominates in the 33—50 keV band. The pulse 
profile changes remarkably in the 50-74 keV band, 
which is the most statistically significant of all the 
energy bands investigated, and is distinguished by a 
broad structure that peaks at around phase 0.0. The 
74—110 keV profile resembles the 50-74 keV one. As 
noted above, the pulse profile above 110 keV is consis- 
tent with random fluctuations. Therefore, our results 
set an observed upper energy bound of 110 keV for 
the hard X-ray pulsations in SGR J1550-5418 during 
this enhanced emission episode. 
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FIG. 4: Evolution of RMS pulsed fraction of 
SGR J1550-5418 as a function of energy. Uncertainties 
are la. The energy bands are the same as those used in 
Figured 



C. Pulsed Fraction 

We computed the RMS pulsed fraction using a 
Fourier based approach as described in 20]. In sum- 
mary, we take the Fourier transform of each pulse pro- 
file, then we calculate the RMS pulsed flux by taking 
the Fourier coefficients of up to third harmonic into 
account, and finally obtain the pulsed fraction val- 
ues by dividing the RMS pulsed flux by the phase- 
averaged flux. In Figure 21 we show the pulsed frac- 
tion spectrum of SGR J1550-5418 in the same energy 
bands as in Figure [31 Although marginally signifi- 
cant, there is an indication of a minimum in the RMS 
pulsed fraction around ^30 keV. The RMS reaches its 
maximum value of 0.55 ± 0.12 in the 50—74 keV band, 
and then dips below detection at energies greater than 
~110 keV. 



IV. SPECTRAL PROPERTIES 

A. Time-Integrated and Time-Resolved Spectral 
Analysis 

We analyzed time-integrated and time-resolved 
spectra of the enhanced emission, using the CSPEC 
Burst data (8.6—897 keV). Similar to the timing anal- 
ysis, we excluded all bursts identified with the untrig- 



gered search within the enhancement period. 

Since the Detector Response Matrices (DRMs) of 
GBM are time dependent due to the continuous slew- 
ing of the spacecraft, a DRM should be generated for 
every 2-3 degrees of slewing (corresponding to every 
~20— 50 s of data). Therefore, we generated DRMs 
for every 50 s starting from To, for this analysis. 

We found clear evidence for spectral curvature be- 
low 100 keV in the time-integrated spectrum of the 
entire burst-free enhancement period (72—248 s); a 
single power law thus resulted in a very poor fit. We 
employed five other spectral models; cut-off power 
law, power law + blackbody, optically-thin thermal 
bremsstrahlung, and single/double blackbody. We 
found that the time-integrated spectrum is best de- 
scribed by a power law + blackbody (see Figure [5]). 
All other spectral models did not provide better fits 
mainly because they failed to fit the lower energy ex- 
cess < 10 keV. The best-fit power-law index and tem- 
perature of the additional blackbody are shown in Ta- 
ble D Adding a blackbody (with kT = 18 ± 4 keV) to 
a power law resulted in the most significant improve- 
ment in Cash statistics [2l[ over a single power law 
(AC-stat = 13.5 for 2 degrees of freedom, correspond- 
ing to an improvement of 3.25ct). 

The average energy flux over the entire enhance- 
ment is (6.5±2.4)xl0~ 8 ergcm- 2 s- 1 (in 8-150 keV), 
of which the blackbody component accounts for 19%. 
We estimate a total isotropic emitted energy of 
2.9 x 10 40 Z? 5 2 erg for the entire persistent emission 
(8—150 keV) during the enhancement. 

To investigate the evolution of the blackbody com- 
ponent and of the source's spectral properties in 
general, we divided the enhanced emission period 
into three time intervals of ~50 s each: 74—117 s, 
122-169 s, and 173-223 s after the trigger time. We 
employed the same set of photon models as the time- 
integrated analysis described above. The first spec- 
trum was best fit by a single power law with no ev- 
idence of a blackbody or any curvature. The second 
and third spectra, on the other hand, were best de- 
scribed by power law + blackbody models. In the sec- 
ond spectrum (the peak of the enhancement) the ad- 
ditional blackbody component was statistically most 
significant, and remained significant in the third spec- 



TABLE I: Spectral parameters of the enhanced persistent 
emission period of SGR J1550-5418. \-a uncertainties are 
shown in parentheses. 



Time 


Power Law 


Blackbody 


Energy Flux 


since To 


Index 


kT 


Power Law Blackbody 


(s) 




(keV) 


(10~ s crgs/cm 2 -s) 


72-248 


-2.06 (0.10) 


17.7 (3.8) 


5.30 (2.37) 1.22 (0.28) 


74-117 


-2.15 (0.17) 


No BB 


2.85 (3.30) 


122-169 


-2.09 (0.11) 


17.4 (1.7) 


7.82 (4.47) 4.08 (0.65) 


173-223 


-2.14 (0.19) 


16.4 (2.7) 


5.05 (3.75) 2.59 (0.72) 



* Flux is calculated in 8—150 keV. 
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FIG. 5: The photon spectrum of the time interval To+122 
to 169 s. The blackbody and power-law components are 
shown separately with dashed curves. The data are binned 
for display purpose only. A 3<j upper limit is shown for 
the highest energy bin. 

trum as well (AC-stat = 42.9 and 15.3, corresponding 
to 6.2(7 and 3.5<r improvements, respectively). The ra- 
tio of the blackbody flux to the total flux (8—150 keV) 
was found to be 34% in both intervals. The indices of 
the underlying power-law component, and the black- 
body temperature also remained constant, at ~— 2.1 
and ~17 keV, respectively (within uncertainties; see 
also Table 1), while the power-law amplitude tracked 
the photon flux. 



B. Phase-Resolved Spectral Analysis 

We performed spin-phase-resolved spectral anal- 
ysis of the pulsed enhanced emission, as follows: 
we co-added the burst-free spectrum of each pulse 
(in To+122 to 223 s, corresponding to second and 
third time-resolved spectra) using TTE data and ex- 
tracted a phase-maximum spectrum (between phases 
0.75—1.25 in Figure [3]) and a phase-minimum spec- 
trum (between phases 0.25—0.75). The spin phase for 
each photon was calculated using barycentered times, 
as was done for the timing analysis. We calculated the 
background spectrum, from the burst-free interval at 
T to T o +60 s. 

The spectra of both the phase minimum and maxi- 
mum were adequately fitted with power law + black- 
body models, where we kept the power-law indices 
and the blackbody temperatures linked. The values of 
the linked parameters found in the fit were consistent 
(within la) with those of the time-integrated spectra 
(see Table HJ. The contributions of the blackbody flux 



to the total flux were (52±18)% and (35±18)%, in 
the phase-maximum and phase-minimum spectra, re- 
spectively. This hints that the blackbody component 
was more significant in the phase-maximum spectrum 
than in the phase-minimum spectrum. 

V. SUMMARY AND DISCUSSION 

We report here the discovery of coherent pulsa- 
tions in the persistent hard X-ray emission from 
SGRJ1550-5418 in the Fermi/ GBM data lasting 
~150 s. Coherent pulsations with a 55% RMS pulse 
fraction have never been detected in the persistent 
emission at these high energies from a magnetar as 
yet. These pulsations were detected only at the on- 
set of a major bursting episode and were not directly 
related to a major burst or flare from the source. Pre- 
viously, intermediate flares with pulsating tails were 
observed from SGR 1900+14 [H[23| and very recently 
from SGRJ1550-5418 (~6 hours after T )fy. Ther- 
mal components were also found in the decaying tails 
of intermediate events from SGR 1900+14 with much 
lower blackbody temperatures of ^2 keV The 
thermal component of the enhanced emission we re- 
port here is hotter (17 keV), exhibits a strong depen- 
dence of the pulse profile with energy with a very high 
RMS pulsed fraction (up to 55%), and is clearly not 
associated with a decaying event tail. Energetically, 
however, the fluence of this enhanced emission is com- 
parable to that of tail emission of the intermediate 
flares from SGR 1900+14. 

Our timing analysis showed that the detection of 
pulsations is most significant in the 120—210 s inter- 
val after trigger. We find that the spectrum requires 
a blackbody component along with a power law be- 
tween 122—223 s, which is consistent with the time 
interval of the most significant detection of pulsa- 
tions. Moreover, as determined by the energy depen- 
dent pulse profiles and RMS pulsed fractions, we find 
that the high-energy pulsations are most significant 
in the 50—74 keV range. Strikingly, the blackbody 
component of the enhanced persistent emission spec- 
trum peaks at around 51 keV (i.e., the Wien peak of 
17 keV, see Figured]). These two independent pieces 
of evidence lend strong support for a blackbody radi- 
ation component to account for the curvature in the 
spectrum of the enhanced emission. 

In our spin-phase-resolved spectral analysis, we 
find that the blackbody flux to the total emission 
is (52+18)% and (35+18)% in the phase-maximum 
and phase-minimum spectra, respectively. This also 
suggests that a major contribution to the observed 
pulsations is from the blackbody component. If we 
assume a surface hot-spot during this pulsating in- 
terval, then the best-fit blackbody corresponds to an 
effective radiating area (as projected on the plane of 
the sky, far from the star) of = ttD 2 F/(<jT 4 ) rj 
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0.046L> 5 2 km 2 , where T « 2 x 10 8 °K (kT w 17 keV) 
is the observed (gravitationally redshifted) tempera- 
ture. We have used here the blackbody flux at the 
peak of the pulsations (i.e., phase maximum; F w 
5 x 10~ 8 erg cm~ 2 s _1 ) where the hot-spot is expected 
to be relatively close to face-on, in order to minimize 
the effects of projection and gravitational lensing by 
the neutron star, so that Soo would be relatively close 
to the physical area, S, of the hot spot on the neutron 
star surface. For a circular hot-spot, this corresponds 
to a radius of ^120D5m. 

The internal field of magnetars can be significantly 
stronger and more tangled than the external dipole 
field. As the internal field twists the stellar crust, 
the magnetosphere also becomes twisted, possibly in 
a complex manner [24| . If we identify the inferred 
size of the hot-spot as the size of a twisted region on 
the stellar surface, the rate of energy dissipation is 
expected to be Ld <C 10 38 erg s _1 for global stabil- 
ity case [25| . This is inconsistent with the observed 
luminosity of the spot, Ld ~ I0 38 L> 5 2 erg s" 1 

It may be possible for the magnetic twist to grow 
to a global instability level during a highly active 
bursting period due to frequent starquakes [25j . As 
the magnetosphere untwists, a large amount of en- 
ergy must be dissipated [2tS| . A small "trapped fire- 
ball" - plasma of pairs and photons confined by 
a closed magnetic field region - could then poten- 
tially account for the inferred hot-spot, and in partic- 
ular its roughly constant temperature and size. Con- 
fining a "fireball" of energy at least comparable to 
that emitted by the observed blackbody component, 
E[ so ,bb ~ 5.6 x 10 39 D 5 2 erg, within a region of radius 
a ~ 120L> 5 m requires Eb(ci) — ^a 3 B 2 > L isoB B or 
B > 1.4 x 10 14 L>~ 1/2 G. This is consistent with the 
surface dipole field ofB~2.2xl0 14 G inferred from 
the measured PP Q. Therefore, a sufficiently small 
closed magnetic loop anchored by the crust could pro- 
vide the required confinement. Moreover, the black- 
body emission is expected to be accompanied by non- 
thermal, high-energy radiation produced by collision- 



less dissipation, with the luminosities comparable to 
the blackbody components [25]. This is in good agree- 
ment with our observations of SGR J1550-5418. 

In conclusion, the area of the blackbody emitting 
region observed here is the smallest "hot spot" mea- 
sured for a magnetar, which likely arises from magnet- 
ically confined hot plasma on the neutron star surface, 
possibly caused by the gradual dissipative untwisting 
of the magnetosphere |26| . If the total radiated en- 
ergy was initially confined to the inferred extremely 
small size of the enhanced emission region (as in a 
mini "trapped fireball" scenario), this would indicate 
a very large magnetic energy density, similar to the 
"trapped fireball" model for the tails of SGR giant 
flares. The observed enhanced emission that we report 
here is much less energetic than a giant flare tail, while 
its energy is comparable to the tail energy of interme- 
diate events and at the high end of typical SGR bursts. 
Despite some distinct properties, the enhanced emis- 
sion of SGR J1550-5418 carries various flavors of all 
three SGR phenomena, and thus it is most likely re- 
lated to the very pronounced bursting activity that 
immediately followed it. 
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